Constraint-Based Qualitative Simulation 
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Abstract 

We consider qualitative simulation involving a finite set 
of qualitative relations in presence of complete knowledge 
about their interrelationship. We show how it can be nat- 
urally captured by means of constraints expressed in tem- 
poral logic and constraint satisfaction problems. The con- 
straints relate at each stage the 'past' of a simulation with 
its future '. The benefit of this approach is that it readily 
leads to an implementation based on constraint technol- 
ogy that can be used to generate simulations and to answer 
queries about them. 



1 Introduction 

Qualitative reasoning was introduced in AI to abstract 
from numeric quantities, such as the precise time of an event 
or the location or trajectory of an object in space, and to 
reason instead on the level of appropriate abstractions. Two 
different forms of qualitative reasoning were studied in the 
literature. The first one is concerned with reasoning about 
continuous change in physical systems, monitoring streams 
of observations and simulating behaviours, to name a few 
applications. The main techniques used are qualitative dif- 
ferential equations, constraint propagation and discrete state 
graphs. For a thorough introduction see 1141 . 

The second form of qualitative aims at reasoning about 
contingencies such as time, space, shape, size, directions, 
through an abstraction of the quantitative information into 
a finite set of qualitative relations. One then relies on com- 
plete knowledge about the interrelationship of these qualita- 
tive relations. This approach is exemplified by temporal rea- 
soning due to Q], spatial reasoning introduced in 1 10 1 and 
1 20 1, reasoning about cardinal directions (such as North, 



°©2005 IEEE. Personal use of this material is permitted. However, 
permission to reprint/republish this material for advertising or promotional 
purposes or for creating new collective works for resale or redistribution 
to servers or lists, or to reuse any copyrighted component of this work in 
other works must be obtained from the IEEE. 



Sebastian Brand 
National University of Singapore, Singapore 
brand @ comp . nu s . edu . sg 



Northwest), see, e.g., 1161 . etc. For a recent overview of 
this approach to spatial reasoning, see JSJ. 

Qualitative simulation deals with reasoning about pos- 
sible evolutions in time of models capturing qualitative in- 
formation. One assumes that time is discrete and that only 
changes adhering to some desired format occur at each 
stage. 1151 discusses qualitative simulation in the first 
framework, while qualitative spatial simulation is consid- 
ered in 1 9 1 . 

Our aim here is to show how qualitative simulation in 
the second approach to qualitative reasoning (exemplified 
by qualitative temporal and spatial reasoning) can be nat- 
urally captured by means of temporal logic and constraint 
satisfaction problems. The resulting framework allows us 
to concisely describe various complex forms of behaviour, 
such as a simulation of a naval navigation problem or a so- 
lution to a version of a piano movers problem. The domain 
knowledge is formulated using a variant of linear temporal 
logic with both past and future temporal operators. Such 
temporal formulas are then translated into constraints. 

The usual constraint-oriented representation of the sec- 
ond approach to qualitative reasoning is based on modelling 
qualitative relations as constraints. See, for example, II II 
for an application of this modelling approach. In contrast, 
we represent qualitative relations as variables. This way 
of modelling has important advantages. In particular, it is 
more declarative since model and solver are kept separate; 
see the study of the relation variable model in (5). In our 
case it allows us to express all domain knowledge on the 
same conceptual level, namely as constraints on the relation 
variables. Standard techniques of constraint programming 
can then be used to generate the simulations and to answer 
queries about them. 

To support this claim, we implemented this approach in 
the generic constraint programming system ECL'PS'' [22 1 
and discuss here several case studies. 



2 Simulation Constraints 

2.1 Constraint Satisfaction Problems 

We begin by briefly introducing Constraint Program- 
ming. Consider a sequence X — x±, . . . , x m of variables 
with respective domains D\, . . . , D m . By a constraint C 
on X, written C(X), we mean a subset of D\ x • • • x D m . 
A constraint satisfaction problem (CSP) consists of a finite 
sequence of variables X with respective domains and a fi- 
nite set C of constraints, each on a subsequence of X. A 
solution to a CSP is an assignment to its variables respect- 
ing their domains and constraints. 

We study here CSPs with finite domains. They can be 
solved by a top-down search interleaved with constraint 
propagation. The top-down search is determined by a 
branching strategy that controls the splitting of a given CSP 
into two or more CSPs, the 'union' of which is equiva- 
lent to (i. e., has the same solutions as) the initial CSP. In 
turn, constraint propagation transforms a given CSP into 
one that is equivalent but simpler. We use here heuristics- 
controlled domain partitioning as the branching strategy 
and hyper-arc consistency of |19| as the constraint prop- 
agation. Hyper-arc consistency is enforced by removing 
from each variable domain the elements not used in a con- 
straint. 

2.2 Intra-state Constraints 

To describe qualitative simulations formally, we define 
first intra-state and inter-state constraints. A qualitative sim- 
ulation corresponds then to a CSP consisting of stages that 
all satisfy the intra-state constraints. Moreover, this CSP 
satisfies the inter-state constraints that link the variables ap- 
pearing in various stages. 

For presentational reasons, we restrict ourselves here to 
binary qualitative relations (e. g., topology, relative size). 
This is no fundamental limitation; our approach extends di- 
rectly to higher-arity relations (e. g., ternary orientation). 

We assume that we have at our disposal 

• a finite set of qualitative relations Q, with a special 
element denoting the relation of an object to itself; 

• consistency conditions on Q-scenarios; we assume the 
usual case that they can be expressed as relations over 
Q, specifically as a binary converse relation conv and 
a ternary composition relation comp, 

• a conceptual neighbourhood relation between the ele- 
ments of Q that describes which atomic changes in the 
qualitative relations are admissible. 



covered by 




covers 



Figure 1. The eight RCC8 relations 

Example. Take the qualitative spatial reasoning with 
topology introduced in [ 1 1 and |20|. The set of qualita- 
tive relations is the set RCC8, i. e., 

Q = {disjoint, meet, equal, covers, 

covered by, contains, inside, overlap}; 

see Fig. [T] which also shows the neighbourhood relation 
between these relations. □ 

We fix now a sequence O of objects of interest. By a 
qualitative array we mean a two-dimensional array Q on 
O x O such that 

• for each pair of objects A, B G O, the expression 
Q[k, B] is a variable denoting the (basic) relation be- 
tween A, B. So its initial domain is a subset of Q. 

• the consistency conditions hold on Q, so for each triple 
of objects A, B, C the following intra-state constraints 
are satisfied: 

reflexivity: Q [A, A] = equal, 
converse: conv(Q[A, B], Q[B, A]), 
composition: comp(Q[A, B], Q[B, C], Q[k, C]). 

Each qualitative array determines a unique CSP. Its vari- 
ables are Q[k, B], with A and B ranging over the sequence of 
the assumed objects O. The domains of these variables are 
appropriate subsets of Q. An instantiation of the variables 
to elements of Q corresponds to a consistent Q-scenario. 

In what follows we represent each stage t of a simulation 
by a CSP Vt uniquely determined by a qualitative array Q t . 
Here t is a variable ranging over the set of natural numbers 
that represents discrete time. Instead of Qt[k, B] we also 
write Q[k, B, t], reflecting that, in fact, we deal with a single 
ternary array. 
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2.3 Inter-state Constraints 

To describe the inter-state constraints, we use as atomic 
formulas statements of the form 

Q[A,B]?g 

where ? G {=,7^} and q G Q, or 'true', and employ a 
temporal logic with four temporal operators, 



□ 



(eventually), 
(from now on), 



O (next time), 
U (until), 



and their 'past' counterparts, O -1 , 0~ , □ , and S 
(since). While it is known that past time operators can be 
eliminated, their use results in more succinct (and in our 
case more intuitive) specifications; see, e.g., 1 18 1. 

Inter-state constraints are formulas that have the form 
— > Oip. Both and -0 are built out of atomic formulas us- 
ing propositional connectives, but contains only past time 
temporal operators and i/j uses only future time operators. 

Intuitively, at each time instance t, each inter-state con- 
straint -> Oip links the 'past' CSP [J* =0 V t with the 'fu- 
ture' CSP Uj=t+i 'Pi ■ ^° we interpret in the interval [0..t], 
and tp in the interval [t + 1 .. i max ]- 

We now explain the meaning of a past or future temporal 
formula with respect to the underlying qualitative array Q 
in an interval [s..t], for which we stipulate s ^ t. We write 
|=r s ..t] to express that holds in the interval [s..t\. 

Propositional connectives. These are defined as ex- 
pected, in particular independently of the 'past' or 'future' 
aspect of the formula. For example, 

\=[s..t] if not \=[ a ..t] 0, 

\=[s..t] 0i V 02 if \=[s..t] (f>i or (=[*..t] <fc- 



Conjunction 0i A 02 and implication 0i 
analogously. 



b% are defined 



Future formulas. Intuitively, the evaluation starts at the 
lower bound of the time interval and moves only forward in 
time. 



\= [a .. t] Q[A,B]?c 

\=[s..t] O0 

\=[s..t] □</> 
\=[s..t] O0 

\=[ s ..t] xU 



if Q[A,B,s]?c 

where ? G {=, 7^}; 

if \=[r..t] <t> 

and r = s + 1 and r ^ t; 
if |=[r..t] for all r G [s..t]; 
if |=[r..t] for some r G [s..i]; 
if Nrr-.t] for some r G [s..i] 
and \= [u .. t ] x for all it G [s .. r - 1]. 
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N[ S ..t] 


Q[A,fl]?c 


if 


N*..t] 


O _1 


if 


N[ S ..t] 


□ -1 


if 


h[ S ..t] 


O _1 


if 


N[*..t] 


XS0 


if 



Figure 2. Navigation path 



Past formulas. Here the evaluation starts at the upper 
bound and moves backward. 



Q[A,B,t}?c 
where ? G {=, 7^}; 

H[s..r] <f> 

and r = t — 1 and s ^ r; 
H[s..r] for all r G [s..t]; 
N[a.-r] for some r G [s..t]; 
H[s..r] for some r G [s..t] 
and (=[„.. t] x for all u G [r + 1 .. t]. 



Furthermore, we write Q[k, B] G {qi, . . . , as an ab- 
breviation of (Q[A, B] = g x ) V ... V (Q[A, B] = q k ). The 
meaning of Q[A, B] ^ {qi, ...,%} is analogous. 

The bounded quantification 3 A G {01, . . . , Ofc}. 0(A) 
represents the disjunction 0(oi) V ... V 0(ofc). Universal 
quantification VA G {01, . . . , o^}. 0(A) is interpreted anal- 
ogously. As usual, A in 0(A) denotes a placeholder (free 
variable), and 0(c>j ) is obtained by replacing A in all its oc- 
currences by Oi. 

2.4 An Example: Navigation 

A ship navigates around three buoys along a specified 
course. The position of the buoys is fixed; see Fig. [2] We 
reason qualitatively about the cardinal directions 

Q= {N,NE,...,W,NW, EQ} 

with the obvious meaning (EQ is the identity relation). 
Ligozat 1 16 1 provides the composition table for this form 
of qualitative reasoning and shows that it captures consis- 
tency. 

The buoy positions are given by the following global 
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intra-state constraints: 

Q[buoy a ,buoy c ] = N\N, 
Q[buoy al buoy b ] = S\N 1 
Q{buoy bl buoy c ] = NW. 

All objects occupy different positions: 

VA, B e O. A ^ B -> Q[A, B] ^ EQ. 

The initial position of the ship is south of buoy c, so we 
have Q[ship, buoy c ] — S. The ship is required to follow a 
path around the buoys. In Fig.|2j the positions required to 
be visited are marked with bold circles. We stipulate 

0(Q[ship,buoy a ] = W A 

0(Q[ship,buoy b } = N A 

<>(Q [ship, buoy c ] = E A 

0(Q[ship,buoy c ] =S)))) 

to hold in the interval [0 .. i max ]. 

A tour of 13 steps exists (and is found by our program); 
it is indicated in Fig. 13 

3 Temporal Formulas as Constraints 

We explain now how a temporal formula (an inter-state 
constraint) is imposed on the sequence of CSPs represent- 
ing the spatial arrays at consecutive times. Such a formula 
is reduced to a sequence of constraints by eliminating the 
temporal operators. We provide two alternative translations. 
The first simply unfolds the temporal operators into primi- 
tive constraints, while the second retains more structure and 
avoids duplication of subformulas by relying on array con- 
straints. 

Consider a temporal formula <p — > Oip where <f> uses only 
'past' time operators and ip uses only 'future' time opera- 
tors. Given a CSP Ui= s we snow now tne P ast temporal 
logic formula <p translates to a constraint cons~ ([s..t], <p) 
and how the future temporal logic formula ip translates 
to a constraint cons + ([s..t], ip), both on the variables of 

u! - a- 

We assume that the target constraint language has 
Boolean constraints and reified versions of simple compari- 
son and arithmetic constraints. Reifying a constraint means 
associating a Boolean variable with it that reflects the truth 
of the constraint. For example, (x — y) = b is a reified 
equality constraint: b is a Boolean variable reflecting the 
truth of the constraint x = y. 

We denote by cons([s..t], <p) = b the sequence of con- 
straints representing the fact that the formula cp has the truth 
value b in the interval [s..t\. The 'past' or 'future' aspect 
of a formula is indicated by a marker ~~ or + , resp., when 
relevant. The translation of cp proceeds by induction and is 
initiated with cons([s..t], cp) = 1 (where s ^ t). 



3.1 Unfolding Translation 

We translate the propositional connectives into appropri- 
ate Boolean constraints. The temporal operators are un- 
folded over the simulation stages. 

For example, the 'future' formula 0(Q[A, B] = q) in the 
interval [1..3] translates to 

(Q[A,B,1] =q) = h, 

(Q[A,B,2]=g) = b 2 , 

(Q[A,B,3] =q) = b 3 , and 
bi V b 2 V 6 3 = 1, 

with fresh Boolean variables b\, 62, 63. 
Translation for 'future' formulas. 



cons + ([s. 


.t], true) = b 


is 


6=1; 




cons + ([s. 


.t],^<P)=b 


is 


V = -6, 










cons + ([s..t], <p) = 




cons + ([s. 


.t],<piV<p 2 ) = b 


is 


(61 V 6 2 ) = b, 










cons + {[s..t], (pi) 


= h 








cons + ([s..t],(p 2 ) 


= h 


cons + ([s. 


.t],Q[A,B]?c) 


= b 


is 






(Q[A,B,s]n 


:) = b where ? 6 {= 


,^}; 


cons + ([s. 


.t],Q<P) = b 


is 


(61 A 6 2 ) = b, 










cons + ([r..t], cp) = 


= 62, 








(s + l^t) = b 1: 










(s + 1 = r) = bi 


1 


cons + ([s. 


.t],D(p) = b 


is 


iKesM = b, 






cons + ([r. 


.t],<f> 


) = b r for all r 6 [s 


••*]; 


cons + ([s. 


.t),o<p) = b 


is 


(V re ,.A) = &. 






cons + ([r. 


.t],4> 


) = b r for all r G [s 


••*]; 


cons + {\s. 


.t], X UcP) = b 


is 








cons + ([s. 


.t],<P 


<VxAO( X U0)) = 


: 6. 



Translation for 'past' formulas. This case is symmetric 
to the 'future' case except for the 'backward' perspective. 
So we have 

cons-([s..t],Q[A,B] ?c) = b is 

{Q[A, B, t] ? c) = b where ? 6 {=, 

for example. The remaining cases are defined analogously. 

Observe that the interval bounds s, t in cons([s..t], cp) are 
treated as constants such that s ^ t. 

3.2 Array Translation 

This alternative translation avoids the potentially large 
disjunctive constraints caused by unfolding the O and U op- 
erators. The idea is to push disjunctive information inside 
variable domains, with the help of array constraints. 
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Reconsider the formula 0(Q[A, B] = q) in the interval 
[1..3]. It is translated into a single array constraint, with the 
help of a fresh variable x ranging over time points: 

Q[k,B,x] = q, 
1 ^ x, x ^ 3. 

Array constraints generalise the better-known element 
constraint. Constraint propagation for array constraints is 
studied in [5 1 and used in our implementation. 

When negation occurs in the formula, a complication 
arises with this translation approach, however. Just negating 
the associated truth value, as in the unfolding translation, is 
now incorrect. We therefore first transform a formula into 
negation normal form (NNF). 

The array translation of NNF formulas follows. We give 
it only for 'future' formulas and where different from the 
unfolding translation. The case of negation does not apply 
anymore. 

cons + ([s..t], 0(j>) = b is 

cons + ([s..t], (j) A (Otrue -> OU<j>) = b; 

cons + ([s..t],O0) = b is s ^ r, r ^ t, 

cons + ([r..t],4>) = b; 

cons+([s..t],x U 4>) = b is (h A {b 2 V b 3 )) = b, 

s ^ r, r ^ t, 
cons + ([r..t],4>) = bi, 
(s = r) = b 2 , 
s u, u ^ r, 
(u = r-l) = b 3 , 
cons + ([s..u], Ox) = b 3 . 

The interval end points s, t in cons([s..t], <f)) can now be 
variables with domains, in contrast to the case of the un- 
folding translation where s, t are constants. We are care- 
ful to maintain the invariant s ^ t and state appropriate 
constraints to this end. Therefore, for example, we un- 
fold \3<f> into a conjunction only step-wise, as the formula 
(j) A (Otrue -> O □</>). 

Example. Let us contrast the two alternative translations 
for a formula from the navigation domain. Consider 

= O(0i A 0^2), 

01 = (Q[ship, buoy] = E) and 2 = (Q[ship, buoy] = S), 
in the interval [l..n] for a constant n, as a 'future' formula. 
So we consider the sequence of constraints cons + ([l..n], <j>) 
for each translation. 

The unfolding translation generates many reified equal- 
ity constraints of the form (Q[ship, buoy, k] = D) = bi t k, 
where D is E or S. More specifically, n + Yn=i * = 
n(n + 3) /2 such constraints and as many new Boolean vari- 
ables are created. Many of the constraints are variants of 
each other differing only in their Boolean variable fc^fc. 



Simulate 

spatial array Q, state constraints, t max 1 — ► solution 

PS := (>;*:= 

while t < t max do 

Vt '■= create CSP from Q t and 

impose intra-state constraints 
VS := append Vt to VS and 

impose inter-state constraints 

{VS, failure) := prop(7 , iS) 
if not failure then 

VS := VS with final state constraints 
imposed on Vt 

(solution, success) := so\ve(VS') 
if success then return solution 
t:=t+l 
return failure 

Figure 3. The simulation algorithm 

The array translation results in just two array constraints, 
namely Q [ship, buoy, ri] — E and Q[ship,buoy,r2\ = S, 
The four ordering constraints 1 ^ r\, r\ ^ n r\ ^ r<z, and 
r 2 ^ n control the fresh variables n, r%. □ 

4 Simulations 

By a qualitative simulation we mean a finite or infinite 
sequence VS = (Vq,Vi, . . .) of CSPs such that for each 
chosen inter-state constraint <j> — > we have that the con- 
straint 

cons ([Q .. to] j <fi) —* cons ( [to + 1 .. t] , ip) 
is satisfied by the CSP |J- =0 V l , 

• if VS is finite with u elements, for alHo £ [0 .. u — 1], 

t t max , 

• if VS is infinite, for all to ^ 0, t ^ to + 1. 

Thus, at each stage of the qualitative simulation, we relate 
its past (and presence) to its future using the chosen inter- 
state constraints. 

Consider an initial situation I = Vq and a final situation 
T x determined by a qualitative array of the form Q x , where 
x is a variable ranging over the set of integers (possible time 
instances). We would like to determine whether a simula- 
tion exists that starts in X and reaches Tt, where t is the 
number of steps. If one exists, we may also be interested in 
computing a shortest one, or in computing all of them. 

Simulation algorithm. The algorithm given in Figure [5] 
provides a solution to the first two problems in presence of 
a non-circularity constraint. 
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The sequence VS of CSPs is initially empty and subse- 
quently step-wise extended; so it remains finite. We view 
VS as a single CSP, which consists of regular finite domain 
variables and constraints and which thus fits into the prob- 
lem format solvable by a standard constraint programming 
techniques. 

We employ the auxiliary procedures prop and solve. The 
call to prop performs constraint propagation of the intra- 
state and inter-state constraints. In our implementation, the 
hyper-arc consistency notion is used. As a result, the vari- 
able domains are pruned and less backtracks arise when 
solve is called. If the outcome is an inconsistent CSP, the 
value false is returned in the failure flag. 

The call to solve checks if a solution to the CSP corre- 
sponding to the given sequence of CSPs exists. If so, a so- 
lution and true is returned, otherwise (0, false). In our im- 
plementation, solve is a standard backtrack search (based on 
variable domain splitting) combined with constraint propa- 
gation as in the prop procedure. 

We use the constant t max equal to the number of different 
qualitative arrays, i.e., i max = \0\ ■ (\0\ - 1) • 21 s !- 1 . If 
the desired simulation exists, the above algorithm finds a 
shortest one and outputs it in the variable solution. 

5 Implementation 

We implemented the simulation algorithm of Fig.|5]and 
both alternative translations of temporal formulas to con- 
straints in the ECL'PS e constraint programming system 
1221 . The total program size is roughly 1500 lines of code. 

5.1 Propagation 

Support for enforcing hyper-arc consistency for Boolean 
and many reified constraints, as well as for extensionally 
defined constraints such as conv, comp and the conceptual 
neighbourhood constraint, is directly available in ECL ! PS e 
(by its fd/ic and propia libraries). For array constraints, we 
use the ECL'PS 6 implementation discussed in |5 1. 

The availability of these (generic) implementations of 
propagation mechanisms explains why we chose hyper-arc 
consistency. We emphasise, however, that in a relation vari- 
able model, constraint propagation is relevant only for effi- 
ciency. 

5.2 Search 

We use the basic backtracking algorithm provided by 
ECL'PS e , but we control it with the heuristics described in 
the following section. 

Various other, advanced search strategies are available 
in ECL'PS e , for example Limited Discrepancy Search 1131 . 
Although we did not experiment with these techniques, we 





H l 
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Figure 4. A piano movers problem 



believe it is worth doing so, and it is not difficult to modify 
our implementation (the solve procedure) accordingly. 

5.3 Heuristics 

Our implementation also incorporates the specialised 
reasoning techniques for RCC8 |21 1 and the cardinal direc- 
tions |16|. In these studies, maximal tractable subclasses 
of the respective calculi are identified, and corresponding 
polynomial decision procedures are discussed. 

Our context requires that these techniques are treated as 
heuristics, due to the presence of side constraints (notably 
the inter-state constraints). With a relation variable model 
for qualitative spatial reasoning, these heuristics fall into the 
customary class of variable and value ordering heuristics for 
guiding search in constraint programming. 

In our implementation, the search heuristic splits the re- 
lation variable domains appropriately so that one of the new 
domains belongs to a maximal tractable subclass of the re- 
spective calculus. 

6 Case Studies 

We now report on two case studies. In both of them, the 
solutions were found by our implementation within a few 
seconds. 

6.1 Piano Movers Problem 

Consider the following version of the piano movers prob- 
lem. There are three rooms, the living room (L), the study 
room (S) and the bedroom (B), and the corridor (C). Inside 
the study room there is a piano (P) and inside the living 
room a table (T); see Figure|4] Move the piano to the living 
room and the table to the study room assuming that none of 
the rooms and the corridor are large enough to contain the 
piano and the table at the same time. Additionally, ensure 
that the piano and the table at no time will touch each other. 
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To formalise this problem, we describe the initial situa- 
tion by means of the following formulas: 



Q[B,L] 


= disjoint A 


Q[B,S] 


= disjoint A 


Q[L,S] 


= disjoint, 


Q[C,B] 


= meet A 


Q[C,L] 


= meet A 


Q[c,s] 


= meet, 


Q[p,s] 


= inside A 


Q[T,L] 


= inside. 



We assume that initially fa, fa, fa hold, i. e., the constraints 
cons~ ([0..0], fa), cons~ ([().. 0], fa) and cons~ ([0..0], fa) 
are present in the initial situation X. 

Below, given a formula <f>, by an invariant built out of <f> 
we mean the formula <f> — » OU<f>. Further, we call a room 
or a corridor a 'space' and abbreviate the subset of objects 
{B, C, L, S} by S. We now stipulate as the inter-state con- 
straints the invariants built out of the following formulas: 

• the relations between the rooms, and between the 
rooms and the corridor, do not change: fa A fa, 

• at no time do the piano and the table fill completely 
any space: 

Vs € S. (Q[P, s] ^ equal A Q[T, s] ^ equal) , 

• together, the piano and the table do not fit into any 
space. More precisely, at each time, at most one of 
these two objects can be within any space: 

Vs G S. -i(Q[P, s] G {inside, coveredby} A 
<2[T, s] G {inside, coveredby}), 

• at no time instance do the piano and the table touch 
each other: 

Q[P,T] = disjoint. 
The final situation is captured by the constraints 

Q[P, L] = inside and Q[T,S] = inside. 

Remarkably, the interaction with our program revealed 
in the first place that our initial formalisation was incom- 
plete. For example, the program also generated solutions 
in which the piano is moved not through the corridor but 
'through the walls', as it were. 

To avoid such solutions we added the following intra- 
state constraints. 

• each space is too small to be 'touched' (met) or 'over- 
lapped' by the piano and the table at the same time: 

Vs G S. ->(Q[a, P] G {overlap, meet} A 
Q[s,T] G {overlap, meet}), 



• if the piano or the table overlaps with one space s, then 
it also overlaps with some other space s', such that s 
and s' touch each other: 

Vs G S. Vo G {P, T}. (Q[s, 0} = overlap -> 

3s' G S. (Q[s', 0} = overlap A Q[s, s'] = meet)), 

• if the piano overlaps with one space, then it does not 
touch any space, and equally the table: 

Vs G S. Vo G {P,T}. 

(Q[s, o] = overlap — > Vs' G S. Q[s', 0} 7^ meet), 

• both the piano and the table can touch at most one 
space at a time: 

Vs,s' G S. Vo G {P,T}. 

(Q[s, 0} = meet A Q[s', 0} = meet — » 

Q[s, s'] = equal). 

After these additions, our program generated the shortest 
solution in the form of a simulation of length 12. In this 
solution the bedroom is used as a temporary storage for the 
table. Interestingly, the table is not moved completely into 
the bedroom: at a certain moment it only overlaps with the 
bedroom. 

6.2 Phagocytosis 

The second example deals with a simulation of phagocy- 
tosis: an amoeba absorbing a food particle. This problem is 
discussed in |9|. We quote: 

"Each amoeba is credited with vacuoles (being 
fluid spaces) containing either enzymes or food 
which the animal has digested. The enzymes 
are used by the amoeba to break down the food 
into nutrient and waste. This is done by routing 
the enzymes to the food vacuole. Upon contact 
the enzyme and food vacuoles fuse together and 
the enzymes merge into the fluid containing the 
food. After breaking down the food into nutri- 
ent and waste, the nutrient is absorbed into the 
amoeba's protoplasm, leaving the waste material 
in the vacuole ready to be expelled. The waste 
vacuole passes to the exterior of the protozoan's 
(i. e., amoeba's) body, which opens up, letting the 
waste material pass out of the amoeba and into its 
environment." 

To fit it into our present framework, we slightly simpli- 
fied the problem representation by not allowing for objects 
to be added or removed dynamically. 
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In this problem, we have six objects, amoeba, nucleus, 
enzyme, vacuole, nutrient and waste. The initial situation 
is described by means of the three following constraints: 

Q [amoeba, nutrient] = disjoint, 
Q[amoeba, waste] = disjoint, 
Q[nutrient, waste] = equal. 

We have the intra-state constraints 

Q[enzyme, amoeba] = inside, 

Q [vacuole, amoeba] 6 {inside, coveredby}, 

Q [vacuole, enzyme] S {disjoint, meet, overlap, covers}, 

and, concerning the nucleus, 

Q[nucleus, vacuole] 6 {disjoint, meet}, 
Q[nucleus , enzyme] 6 {disjoint, meet}, 
Q[nucleus, amoeba] = inside. 

The inter-state constraints are 

Q[nutrient, amoeba] = meet — > 

O Q[nutrient, amoeba] = overlap, 

Q[nutr., amoeba] 6 {inside, coveredby, overlap — > 

O Q[nutr. , amoeba] G {inside, coveredby}. 

We model the splitting up of the food into nutrient and waste 
material by 

Q[nutrient, waste] = equal -A 

(</>! 4> 2 v 4> 3 ) 
v 

O Q[nutrient, waste] ^ equal; 

with 

4>i = Q[nutrient, vacuole] = inside A 
Q [enzyme, nutrient] = overlap A 
Q[enzyme, waste] = overlap 

<f>2 = O Q[nutrient, waste] = overlap 

03 = O Q[nutrient, waste] = equal 

The dotted operators express if-then-else, that is, 

a -A bV c = (a — > 6) A (~<a — » c). 

The final situation is described by means of the con- 
straints 

Q [amoeba, waste] = disjoint, 
Q[amoeba, nutrient] 6 {contains, covers}. 

Our program generated a simulation consisting of 9 steps. 



7 Final Remarks 

The most common approach to qualitative simulation is 
the one discussed in 1 14 chapter 5]. For a recent overview 
see 1151 . It is based on a qualitative differential equation 
model (QDE) in which one abstracts from the usual differ- 
ential equations by reasoning about a finite set of symbolic 
values (called landmark values). The resulting algorithm, 
called QSIM, constructs the tree of possible evolutions by 
repeatedly constructing the successor states. During this 
process, CSPs are generated and solved. 

This approach is best suited to simulate evolution of 
physical systems. A standard example is a simulation of 
the behaviour of a bath tub with an open drain and constant 
input flow. The resulting constraints are usually equations 
between the relevant variables and lend themselves natu- 
rally to a formalisation using CLP(FD), see (7J chapter 20] 
and |3|. The limited expressiveness of this approach was 
overcome in |4|, where branching time temporal logic was 
used to describe the relevant constraints on the possible evo- 
lutions (called 'trajectories' there). This leads to a modified 
version of the QSIM algorithm in which model checking is 
repeatedly used. 

Our approach is inspired by the qualitative spatial sim- 
ulation studied in |9|, the main features of which are cap- 
tured by the composition table and the neighbourhood re- 
lation discussed in Example 12.21 The distinction between 
the intra-state and inter-state constraints is introduced there, 
however the latter only link the consecutive states in the 
simulation. The simulation algorithm of |9| generates a 
complete tree of all 'evolutions', usually called an envision- 
ment. 

In contrast to |9|, our approach is constraint-based. This 
allows us to repeatedly use constraint propagation to prune 
the search space in the simulation algorithm. Further, by us- 
ing more complex inter-state constraints, defined by means 
of temporal logic, we can express substantially more so- 
phisticated forms of behaviour. 

While the prevalent approach to constraint-based mod- 
elling of qualitative spatial knowledge maps qualitative re- 
lations to constraints, we use variables to express qualita- 
tive relations. The relation variable approach is much more 
declarative, separating the model from the solver. The ad- 
vantage of a relation variable model for qualitative simu- 
lations is that the knowledge of the spatial domain as well 
as of the application domain can be expressed on the same 
conceptual level, by intra-state and inter-state constraints. 
This leads to a model that can easily be realised within a 
typical constraint programming system using generic prop- 
agation and search techniques, and is also immediately open 
to advances in these systems. 

Simulation in our approach subsumes a form of plan- 
ning. In this context, we mention the related work 1 17 1 in 
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the area of planning which shows the benefits of encoding 
planning problems as CSPs and the potential with respect 
to solving efficiency. Also related is the TLPLAN system 
where planning domain knowledge is described in tempo- 
ral logic |2|. The planning system is based on incremental 
forward-search, so temporal formulas are just unfolded one 
step at a time, in contrast to the translation into constraints 
in our constraint-based system. 

Finally, |12| discusses how a qualitative version of the 
piano movers problem can be solved using an approach 
to qualitative reasoning based on topological inference and 
graph-theoretic algorithms. Our approach is substantially 
simpler in that it does not rely on any results on topology 
apart of a justification of the composition table. 
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